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Abstract: 

General  forms  for  the  eigenvalues  and  eigenvectors  of  certain  patterned 
correlation  matrices  are  obtained-  The  results  are  discussed  in  the  context  of 
a  factor  analysis  of  common  stock  rates  of  return. 
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1.  Introduction 

Me  derive  the  general  forms  for  the  eigenvalues  and  eigenvectors  of  certain 
correlation  matrices  with  special  structures.  Our  inquiry  was  motivated  by  an 
analysis  of  the  rates  of  return  on  common  stocks.  Stocks  within  certain  groups 
(for  example,  stocks  within  the  same  industry)  may  be  equally  correlated  and  the 
correlations  between  stocks  in  different  groups,  although  equal,  may  be  smaller 
in  magnitude  than  the  within  group  correlations. 

The  correlation  matrices  we  examine  are  composed  of  blocks  of  equally  cor¬ 
related  groups.  The  correlations  within  each  block  are  expressed  as  simple 
powers  (0,1,2  and  3)  of  a  basic  correlation  p.  For  example,  with 


we  can  form  the  "two  group"  correlation  matrix: 

P 

(m+n)x(m+n) 

where  and  A2  are  analogous  to  A  with  dimension  (mxm)  and  (nxn)  respectively. 

In  the  following  sections  we  consider  several  patterned  correlation 
matrices  of  the  type  discussed  above.  Our  eigenvalue-eigenvector  results  can  be 
generally  applied  to  many  problems  involving  a  principal  component  or  factor 
analysis.  However,  it  will  be  convenient  to  refer  to  the  groups  of  variables  as 
groups  of  common  stock  rateu  of  return  and  to  interpret  our  work  in  terms  of  the 


factor  structure  for  stock  prices  embedded  in  the  Arbitrage  Pricing  Theory  (APT) 
of  assets  proposed  by  Ross  (1976).  We  provide  a  small  example  of  factor  struc¬ 
ture  with  five  stocks  from  two  industries. 


Specifically,  we  consider  correlation  matrices  corresponding  to  the  follow¬ 
ing  cases: 

1.  Two  groups  of  equally  correlated  variables  (stocks)  with  m=n  (equal 
group  sizes)  and  m*n  (unequal  group  sizes) 

2.  Three  groups  of  equally  correlated  variables  (stocks)  with  m=n=p  (equal 
group  sizes) 

3.  Four  groups  of  equally  correlated  variables  with  m=n=p=q  (equal  group 
sizes) 


2. 


Patterned  Correlation  Matrices  for  Two  Equally  Correlated  Groups 
We  define 

1  p  ...  i 


A  = 


P  1 


P  P 


1 


B  = 


2  2 
P  P 


i1 


2  2 
P  P 


n2  2 
P  P 


(1 


with  -1  <_  p  <_  1.  Our  analysis  is  restricted  to  p  >  0  and  all  of  our  examples 
deal  with  this  case.  The  matrix  A  is  square  with  dimensions  determined  by  the 
problem.  In  this  section  the  dimensions  of  A  are  (mxm)  or  (nxn).  The  matrix  B 
is,  in  general,  rectangular.  Here  B  has  dimension  (mxn). 


Groups  of  equal  sizes:  m=n 

Using  (1),  consider  the  (2mx2m)  correlation  matrix 
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The  matrix  p  implies  that  stocks  within  the  two  individual  groups  are  equally 

correlated  with  correlation  coefficient  p.  The  correlation  between  any  two 

2 

stocks  from  different  groups  is  p  . 

Making  use  of  the  chosen  pattern  for  p,  it  can  be  shown  that  p  has  three 
distinct  eigenvalues: 

1 ^  =  1  -  p  (multiplicity  2(m-1)) 

(2J  2 

A1  #  «  1  +  (m-1 )  p  -  mp  (3) 

X(3)  =  1  +  (m-1)  p  +  mp2 

with  X(3)  >  A^2^  >  A* ^  ^  0  (since  by  assumption  p  >  0)*  The  2(m-1)  eigenvec 

tors  (unnormalized)  corresponding  to  the  smallest  eigenvalue  A^1^  =  1  -  p  are 
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The  eigenvectors  (unnormalized)  corresponding  to  the  remaining  eigenvalues  are: 
For  *(2)  =  1  +  (m-1)p  -  rap2. 


V2ra-1  "  I1'  1f  “*'  1f  -1'  _1'  *•"  "11’ 


(5) 


(2)  2 
For  A1  =  1  +  (m-1)p  +  mp  , 


[l,  1,  •  •  •  i  1  ] * 


If  we  regard  the  eigenvectors  as  (unsealed)  coefficient  vectors  for  the  princi¬ 
pal  components  obtained  from  the  correlation  matrix  p,  then  we  see  that  the 
eigenvectors  corresponding  to  the  two  largest  eigenvalues  represent  a  "sum"  of 
all  the  stocks  (variables)  and  a  "contrast"  between  the  stocks  (variables)  in 
the  two  groups  respectively.  A  similar  interpretation  of  the  eigenvectors  can 
be  given  in  terms  of  an  underlying  factor  model  (see  section  5). 


Group  of  unecfual  sizes:  m*n 

Consider  the  correlation  matrix  p  in  (2)  with  B  of  dimension  mxn,  m*n  and 

“w 

the  second  diagonal  A  matrix  of  dimension  (nxn).  It  is  straightforward  to  verify 
that  one  of  the  distinct  eigenvalues  of  p  is  =  1  -  p  with  multiplicity 

(m-1)  +  (n-1)  =  m+n-2.  The  m  +  n-2  (unnormalized)  eigenvectors  corresponding 
to  A  =1-0  are; 


[1*  1  #  1#  •••# 


—  ( m- 1)?  0,  •  •  •  f 


0J’ 


(6) 


1 


V 

m 


[o,  0?  1/  “1  /  0,  •  •  •  i  0 j* 


VnH-n-2~  t0'  *•*'  0;  1f  1#  1#  •**' 

Tbe  remaining  two  eigenvectors  and  associated  distinct  eigenvalues  can  be 
determined  from  orthogonality  requirements.  Specifically,  the  remaining  eigen¬ 
vectors  must  be  of  the  form  V  =  [a,  a,  ...»  a;  b,  b,  b]'  with  both  a*0  and 

b*0.  From  the  definition  of  eigenvalues  and  eigenvectors  we  thus  arrive  at  two 
distinct  equations: 

a(1  +  (m-1)p)  +  b(np2)  =  X  a 

a(mp2)  +  b(  1  +  (n-1)p)  =  Xb 

We  set  a  =  1  without  loss  of  generality  so  that 

2 

1  +  (m-1 ) p  +  bnp  =  X 
5^-  +  1  +  (n-1 )  p  =  X 
or,  equivalently. 


2 

1  +  (ra-1)p  +  bnp  =  1  +  (n-1)p  + 


2 


(7) 


.•jej 


b  = 


+  /  2  2~ 
-(m-n)p  —  p  /(m-n)  +  4mnp 


—  Xm-n) 


2np 


-(m-n)  —  /(m-n)  +  4mnp 

2np 


(8) 


Denoting  the  first  (largest)  solution  for  b  in  (8)  by  b+  and  the  second 
(smallest)  solution  for  b  by  b”  we  obtain  the  remaining  two  distinct  eigenvalues 
and  corresponding  (unnormalized)  eigenvectors.  Thus 

(2)  2 
X'  =  1  +  (m-1)p  +  b  (np  ) 

Vm+n-1  =  1#  *'  11  b  '  b  .  b  J* 

X(3)  =  1  +  ( m- 1 )  p  +  b+  (np2)  (9) 

Vn  =  f1'  1f  1f  b+'  b  '  b 

Noting  that  b  b  - - ,  it  is  clear  that  V  ,  ,  and  V  are  orthogonal  and 

n  ra+n- 1  m+n 

each  of  these  vectors  is  orthogonal  to  the  eigenvectors  vm+n_2*  tore- 

<v  *v 

over,  since  b+  >  0  and  b+  >  b  ,  X^3^  >  X^.  It  can  also  be  shown  that 

x(2)  -  X(1)  >  0  and  this  difference  increases  as  m  and/or  n  increases.  Thus  we 

have  X(3)  >  X^2^  >  X^1^  and  as  m  and/or  n  increases  the  "importance"  of  the 

principal  components  constructed  from  the  eigenvectors  corresponding  to 

.(3)  „  .(2)  , 

X  and  X  increases. 

We  note  that  the  elements  of  V  ,  ,  the  eigenvector  associated  with  the 

m+n 

largest  eigenvalue,  are  all  positive  and  become  equal  as  p  +  1.  The  first  m 

elements  of  the  eigenvector  vm+n  ate  of  different  sign  than  the  last  n 

elements.  In  addition,  the  sum  of  the  elements  of  V  .  ,  approaches  zero  as 

m+n- 1 

p  +  1.  Again,  from  a  principal  components  viewpoint,  if  p  is  positive  and 
large,  the  eigenvectors  corresponding  to  the  two  largest  eigenvalues  represent  a 


sum"  and  a  "contrast"  of  the  original  m+n  variables 


Let  m=3,  n=2  and  p=. 6  then 


A1 

(3x3) 

B 

(3x2) 

1.00 

0.60 

0.60 

0.60 

1.00 

0.60 

0.60 

0.60 

1.00 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 

1.00 

0.60 

B* 

(2x3) 

A2 

(2x2) 

0.36 

0.36 

0.36 

0.60 

1.00 

From  (8), 


-1  +  /l  +  4(6) (0.367 

2(2) (0.6) 


0.8770 


b 


-1  -  /l  +  4(6) ( 0. 36T 

2(2) (0.6) 


-1.7103 


Thus  X*1*  =  1  -  p  =  0.4  (multiplicity  3),  and  X^,  X^  can  be  determined  using 
b+  «  0.8770,  b"  »  -1.7103  and  (9).  We  have: 


Eigenvalues  (Ordered  according  to  decreasing  magnitudes) 


2.8315 

0.9685  0.4000 

0.4000 

0.4000 

Normalized 

eigenvectors  (columns) 

0.4694 

0.3361 

0.7071 

0.4082 

0.0000 

0.4694 

0.3361 

-0.7071 

0.4082 

0.0000 

0.4694 

0.3361 

0.0000 

-0.8165 

0.0000 

0.4117 

-0.5749 

o.oooo' 

0.0000 

0.7071 

0.4117 

-0.5749 

0.0000 

0.0000 

-0.7071 

3.  Correlation  Matrix  for  Three  Equally  Correlated  Groups 

We  consider  three  groups  of  stocks  of  the  same  sizes,  so  m=n=p.  Using  the 


definitions  in  (1),  consider  the  (3mx3m)  correlation  matrix 
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where  both  A  and  B  are  of  dimensions  (mxm).  Thus  stocks  within  each  of  the 

three  groups  are  equally  correlated  with  correlation  coefficient  p.  The  corre- 

2 

lation  between  any  two  stocks  with  each  stock  from  a  different  group  is  p  . 

For  the  matrix  p  in  (10)  there  are  three  distinct  eigenvalues.  The  smallest 
eigenvalue  A^  =  1  -  p  has  multiplicity  3(m-1)  and  the  associated  (unnormal¬ 
ized)  eigenvectors  can  be  written  as 


V.j  —  [lr  -1/  Of  »..f  0;  Of  ...f  0;  Of  ...f  oj* 

V0  =  [  1 ,  If  -2f  ...f  0;  0,  « .  •  f  0;  0 ,  0  ]  * 


—  [if  1#  1#  . • • f  — ( m— 1 )  j  Of  ...f  0#  Of  • • • f  0] 

n—  I 

=  [Of  ...f  0;  1,  -1,  0,  ...f  0;  0,  •  «.,  0  ] ' 


V  =  [o,  ...f  0;  If  If  If  »..f  — (m— 1);  Of  •••f  0 ] * 


V2m-1  "  t° . °!  ° . °?  U  °'  ***' 


V3m-3  =  t° . °;  ° . °?  h  1 . '(m‘1)]’ 

Using  the  orthogonality  requirement,  the  remaining  three  eigenvectors  have 
the  form  V  =  [a,  a;  b,  ...»  b;  c,  c].  Assume  first  that  a*0  (the  case 

of  a=0  will  be  discussed  later)  and  hence,  without  loss  of  generality,  take  a=1. 
Solving  p  V  =  AV  with  V  =  [l,  . 1;  b,  b;  c,  c]  yields  three 

separate  equations 


1  +  (m-1)p  +  mp  b  +  mp  c  =  X 
«  2 

mp  +  (1  +  (m-1)p)b  +  mp  c  =  Xb  (12) 

2  2 

mp  +mpb+(1  +  (m-1)p)c  =  Xc 

Multiplying  the  first  equation  in  (12)  by  b  and  subtracting  the  second  equation 
gives 

(b2  -  1)mp2  +  mp2c  (b  —  1 )  =  0 

A  similar  manipulation  using  the  first  and  third  equations  in  (12)  gives 

(c2  -  1)mp2  +  mp2b(c  -  1)  =  0 
2 

Since  mp  *  0  (for  p*0) ,  we  can  write 
b2  -  1  +  c(b  -  1)  =  0 

c2  -  1  +  b(c  -  1)  =  0  (13) 

or,  after  subtracting  the  second  equation  immediately  above  from  the  first, 
b2  -  c2  +  (b  -  c)  =  0 

Consequently,  (b  -  c) (b  +  c  +  1)  =0  and  hence  b  =  c  or  b  +  c  +  1  —  0. 

Case  I:  Suppose  b  =  c.  Solving  the  first  equation  in  (13)  gives 

b  =  c  =  1  or  b  =  c  =  -1/2 

Case  XI:  Suppose  b  +  c  +  1  =  0.  Any  combination  of  b  and  c  satisfying 

b  +  c  +  1  =  0  solves  (13)  and  consequently  there  are  an  infinite 
number  of  solutions. 

To  summarize,  with  a  =  1,  the  condition  b  =  c  =  1  gives  the  largest  eigen— 
(3) 

value,  X  ,  and  corresponding  eigenvectqr 
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f  31  2 

X  =  1  +  (m  -  1)p  +  2m  p 


V3m  =  f1'  **•'  1f  ^  •**'  1;  1'  **•'  1]* 


(14) 


Using  (12)  with  b  =  c  =  -1/2  or  b  +  c  +  1  =  0,  we  get  the  second  largest  eigen- 
( 2 ) 

value,  X  “  ,  and  two  associated  eigenvectors 


,(2)  ,  A  ,  ,,  2 

X  =  1+  (m-  1)p-mp 


^3m~2  =  [^'  •••/  If  “1/2,  •  •«,  —1/2;  —1/2,  •••»  — 1/2] 


(15) 


V3m-1  t1'  •**'  1f  b'  •**'  b;  c'  **•'  cl’ 

«%<« 

with  b  +  c  +  1  =  0 

Note  the  vectors  in  (15)  are  not  orthogonal. 

Starting  with  b  =  1  or  c  =  1  gives  the  following: 

b  =  1 


(3)  2 

X  =1+(m-1)p+  2m p 


3m 


£  1 ,  •  • . ,  1;  1,  •  •  •  ,  1;  1,  ...,  1  j * 


X(2>  =  1  + 


1  +  (m-1 )  p  -  mp 


(16) 


V3m-2  =  f-1/2'  **•'  -1/2;  1#  •••/  If  “1/2,  ...»  -1/2]' 


V3m-1  "  ta'  ***'  a?  ***'  17  c'  ***'  c^’ 


with  a  +  c  +  1  =  0 


A'  =1  +  (m-1)p  +  2mp 


V  [l,  ...»  1;  1,  1;  1,  1  ] ' 


A^2^  =  1  +  ( m- 1 ) p  -  mp2 


—  [~1/2/  •••/  —1/2;  —1/2/  ♦••/  —1/2;  1/  •*•/  ij1 

V3m-1  =  ^a'  a;  b'  b;  1f  **•'  11’ 


with  a  +  b  +  1=0 


Combining  the  results  in  (15),  (16)  and  (17)  we  see  that  some  of  the  eigen- 

(2)  2 

vectors  corresponding  to  the  eigenvalue  A  =1+(m-1)p-mp  can  be  written 
as  linear  combinations  of  the  other  eigenvectors  corresponding  to  this  eigen¬ 
value.  If  we  exclude  the  redundant  eigenvectors  and  orthogonalize  the  remaining 
two  eigenvectors  we  have  the  final  solution: 

(3)  2 

A  =  1  +  (m-1 ) p  +  2mp 


V3m  =  f1'  **•'  1f  1r  ***'  1;  1f  1J’ 


(2)  2 
A  -  1  +  (m-1)p  -  mp  (multiplicity  2) 


V3m-2  =  f3/2'  *•*'  3/2j  °'  •**'  0;  "3/2,  —3/2 ]* 

V3m-1  =  H/2'  •**'  "V2;  1#  1;  "1/2f  •••'  _1/2]’ 


Recall  our  discussion  leading  to  (18)  started  with  a=1  (and  b=1  and  c=1). 

It  is  easy  to  see  that  an  initial  choice  a=0  (or  b=Q  or  c=0)  leads  directly  to  a 


particular  case  already  considered  for  b=1  (or  c=1  or  a=1).  Therefore,  we  do 


12 


not  have  to  investigate  the  choice  a=0  (or  b=0  or  c=0)  and  (18)  does  indeed 
represent  the  two  largest  distinct  eigenvalues  and  their  eigenvectors. 


A  numerical  example  (principal  components  and  factor  analysis 

Let  m=3  and  p=.7  so  the  correlation  matrix  (10)  is  (9x9).  Using  (18)  and 
A(1)  =  1  -  p,  the  eigenvalues  of  p  are: 


A.  *  1  +  2p  +  6p  =  5.34 


A2  =  Aj  =  1  +  2p  -  3p  = 


=  .93 


A  —  A  ...  —  A  —  .3 
4  5  9 


Let  Z  denote  the  (9x1)  vector  of  standardized  variables. 

Let  e  ,  e„,  . ..,  eQ  denote  the  unit  length  eigenvectors  corresponding  to  the 

1  c.  9 

*v<  V 

eigenvalues  A^,  ...,  Ag.  The  first  principal  component,  Y^  =  e.j'Z,  explains 

-w  <v 

a  proportion  A^trace  ( p)  =  5.34/9  =  .539  of  the  total  variance  of  Z  (see  John- 

*N# 

son  and  Wichern,  1982,  chapter  8).  The  cumulative  proportion  of  total  variance 
explained  by  the  first  two  principal  components,  Y.^  =  e^'Z  and  Y2  =  e2'Z,  is 

*V  -V  <v 

(5.34  +  .93)/9  =  .697  and  so  on. 

Consider  the  factor  model  (see  Johnson  and  Wichern,  1982,  chapter  9) 


Z  =  L  F  +  e 


(9x1)  (9xk)  (kxl)  (9x1) 


where  E(F)  =  0;  Cov(F)  =  I;  E(e)  =  0;  Cov(e)  =  ¥  (diagonal);  E(cF')  =  0  and  L 
is  a  matrix  of  factor  loadings.  Let  k=1.  Representing  the  loadings  by  the 
scaled  eigenvector  e1  (see  Johnson  and  Wichern,  1982),  we  have  L'  =  /A^  e^'  = 

^  *v 

[.770,  .770,  ...,  .770  ].  The  communalities  are  h^  =  .593,  i  =  1,  2,  ...»  9  and 
hence  the  specific  variances  are  ¥  =  .407,  i  =  1,  2,  ...,  9.  The  factor  model 


gives 


Cov(Z)  =  I,L'  +  S’  = 


(symmetric) 


and  the  sum  of  squared  errors  in  the  "residual  matrix,"  p  -  (LL*  +  ¥),  is  .779.  Hie 
performance  of  k=2  and  k=3  factor  models  can  be  analyzed  in  a  similar  fashion.  Note 
for  a  k=2  factor  model  there  is  a  choice  for  the  loadings  on  the  second  factor  since 
^2  ~  and  consequently,  L  =  [t^e^,  /X^e2]  or  L  =  ^3e3]. 

4.  Correlation  Matrix  for  Four  Equally  Correlated  Groups 

We  consider  four  groups  of  stocks  of  the  same  sizes  so  m=n=p=q.  Using  the 
definitions  in  (1)  and  setting 


C  = 

( mxm)  .  .  .  . 


we  consider  the  (4mx4m)  correlation  matrix 


A_ 

B 

C 

C 

B_ 

A 

C 

C 

C 

C 

A 

B 

C 

C 

B 

A_ 

where  A,  B  and  C  are  all  of  dimension  (mxm).  Thus  we  have  identified  two  sets 
of  groups  (1,2),  say,  and  (3,4). 

If  we  regard  the  groups  as  stocks  organized,  for  example,  according  to 
industries,  then  stocks  within  the  same  industries  are  equally  correlated  with 


the  same  correlation  p  for  all  industries.  Stocks  for  industries  1  and  2 
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respectively,  are  equally  correlated  with  correlation  p  .  A  similar  result 
holds  for  industries  3  and  4.  Stocks  in  different  industry  groups,  (1,2)  and 
(3,4),  respectively,  are  equally  correlated  with  correlation  p3. 

Given  the  structure  of  p  in  (21),  it  is  easy  to  see  that  one  of  its  eigen¬ 
values  is  =  1  -  p  with  multiplicity  4(m-1).  The  eigenvectors  are  of  the 

same  form  as  (4),  (6)  and  (11)  so  we  have 


=  1  -  p 


C -  s*  /  ■»  /  '  - «■  ■«  "  o 

=  [1*  -1,  •••,  0;  0,  0;  0,  .»•,  0;  0,  ...,  o]1 


—  [ 1 ,  1,  ~ ( m— 1);  0,  •  •  • ,  0;  0,  •••,  0;  0,  • • » ,  0  ] 1 

=  [0,  0;  1,  -1,  ...,  0;  0,  ...»  0;  0,  0 ] ' 


V2(m-1)  "  I0'  •**'  0?  1f  1f  **•'  “(TO_1);  °»  0;  °»  °I * 

V2m-1  =  f°'  ***'  0;  °'  **•'  0f  1#  **1#  ***'  0f  °'  •**'  °]' 


3(m-1) 


c  •••/  0;  0 f  •••/  0/  1/  1/  •••/  •*( m- 1);  0/  •••/  0  ] ' 


[0/  •••/  0/  0/  •••/  0/  0/  •••/  0/  1/  *1/  • • • /  0 j* 


4(m-1) 


=  [0,  •••#  0;  0,  •••#  0;  0,  •  ••/  0;  1,  •  ••#  -(m-1)]' 


The  remaining  four  eigenvectors  must  be  of  the  form 


V—  [a*  •••#  ei?  b,  •••#  b?  c /  c?  d/ 

equivalently , 

2  3  3 

(1  +  (m-1)p)a+mp  b  +  mp  c  +  mp  d  =  Xa 

2  3  3 

mp  a  +  (1  +  (m-1)p)b  +  mp  c  +  mp  d  =  Xb 

3  3  2 

mp  a  +  mp  b  +  (1  +  (m-1)p)c  +  mp  d  =  Xc 

3  3  2 

mp  a  +  mp  b  +  mp  c  +  (1  +  (m-1)p)d  =  Xd 


and  must  satisfy  p  V  =  XV  or 


(23) 


Algebraic  manipulations  of  the  equations  in  (23),  along  the  lines  of  the  manipu¬ 
lations  in  section  3,  yield  the  remaining  distinct  eigenvalues  and  associated, 
orthogonal,  eigenvectors.  These  eigenvalues  and  (unnormalized)  eigenvectors  ares 


*(4) 


2  3 

=  1  +  ( m- 1 ) p  +  mp  +  2mp 


x(3) 

V4m-1 

*(2) 


V 

4m-2 

v 

4m-3 


—  |l,  ...,  1 j  1,  ...,  1 1  1,  ...,  1?  1,  •#.,  1  ] 

2  3 

=  1  +  (m-1)p  +  mp  -  2mp 

=  [l,  •••,  1;  1,  •••,  1;  ~1,  —l]' 

2 

=  1  +  (m-1)p  -  mp  (multiplicity  2) 

3  •  •  • ,  1 1  —  1 ,  •  •  • ,  “h  0,  •••,  0?  0,  •••,  0  j ! 

=  [o,  ...»  0;  0,  ...»  0;  1,  •••,  1;  -1,  -l]' 


(24) 


Regarding  the  eigenvectors  as  the  coefficients  in  principal  components  con¬ 
structed  from  standardized  variables,  we  "see  that  the  principal  component  corre- 

(4) 

sponding  to  the  largest  eigenvalue,  X  ,  is  a  "sum"  of  the  variables;  the  prin¬ 
cipal  component  corresponding  to  X^3^  represents  a  "contrast"  between  the  groups 
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(2 ) 

(1*2)  and  (3,4)  and  the  principal  components  corresponding  to  X  contrast  col¬ 
lections  1  and  2  and  collections  3  and  4  respectively. 

Suppose  the  sample  correlation  matrix  for  the  rates  of  return  on  m  common 
stocks  from  each  of*  say*  four  different  industries  is  of  the  form  (21).  Assume  the 
factor  model  (19)  with  k=4.  Taking  the  loadings  matrix  L  to  be  L  =  [vT.e  , 

^e4^  Where  V  X<4)'  el  =  V4m/,V4m1'-  X2  =  e2  =  V4m-/"V4m-1  “5 

V  >  V  V  V  V  V 

( 21 

X  =  X.  =  A'  ,  e  =  V  /8V  1  and  e,  =  V(  „/lV  I;  we  might  identify  the 

3  4  3  4m- 2  4m- 2  4  4m-3  4m-3 

*W  *N*  «V  V  N.  V 

first  factor  as  a  general  market  factor,  the  second  factor  as  a  bipolar  factor 
contrasting  industries  (1,2)  with  industries  (3,4).  The  third  factor  as  a  fac¬ 
tor  contrasting  industries  1  and  2  and  the  fourth  factor  as  a  factor  contrasting 
industries  3  and  4.  Interpretations  of  this  kind  for  stock  prices  have  been 
suggested  by  King  (1966).  These  interpretations  are  also  consistent  with  the 
Arbitrage  Pricing  Theory  (APT)  of  assets  suggested  by  Ross  (1976). 

5.  Patterned  Correlation  Matrices  and  Stock  Price  Behavior 

The  APT  of  Ross  (1976)  starts  with  an  assumption  on  the  stock  return 
generating  process.  The  return  on  the  ith  stock,  X^,  is  assumed  to  be  generated 
by  a  k  factor  linear  model  of  the  form: 


X. 

x 


U. 

x 


*  *ii  p< + 


+  F  +  e. 

ik  k  x 


i  =  1,  2,  ...,  N 


(25) 


The  first  term,  ix,  is  the  expected  return  on  the  ith  stock.  The  k  factors, 

F.j,  are  common  to  all  stocks  under  consideration.  The  coefficients  JL.. 

(factor  loadings)  measure  the  reaction  of  stock  i's  returns  to  movements  in  the 
common  factors  Fj .  The  common  factors  represent  the  systematic  risk  components 
in  the  model.  The  error  term,  e^,  represents  the  unsystematic  risk  component 


specific  to  the  ith  stock.  Finally,  it  is  assumed  that  N,  the  number  of  stocks 


(assets),  is  much  greater  than  k,  the  number  of  underlying  common  factors.  If 
the  expected  returns  are  subtracted  from  the  returns  and  the  risk  components  are 
rescaled,  the  model  in  (25)  is  analogous  to  the  linear  factor  model  (19). 

Beginning  with  the  linear  factor  model  for  returns,  Ross  develops  an  equi¬ 
librium  "pricing"  relationship  for  stocks  (assets)  which  must  hold  in  the 
absence  of  riskless  arbitrage  opportunities.  The  accompanying  APT  has  achieved 
a  significant  level  of  acceptance  in  the  finance  academic  community.  It  has 
been  tested,  empirically,  with  varying  degrees  of  support.  In  particular.  Roll 
and  Ross  (1980)  claim  to  find  support  for  the  APT  in  their  empirical  investiga¬ 
tion  and  suggest  the  number  of  common  factors  is  around  four  or  five.  It  is  the 
latter  result,  along  with  our  frequently  observed  "regular"  correlations  between 
stock  returns,  that  motivated  our  investigations  in  the  previous  sections. 

As  a  small  "real  world"  example,  Johnson  and  Wichern  (1982,  chapter  8)  report 
the  correlation  matrix,  R,  constructed  from  one  hundred  daily  rate  of  returns  for 
each  of  five  stocks.  Three  of  the  stocks  are  chemicals  (Allied  Chemical,  DuPont  and 
Union  Carbide)  and  two  of  the  stocks  are  oils  (Exxon  and  Texaco).  The  correlation 
matrix  and  its  eigenvalues  and  normalized  eigenvectors  are  given  below. 

chemicals 

~ — - >  — 

1.0  .577  .509  .387  .462 

1.0  .599  .389  .322 

_ 1U0 _ .436  .426 

/  „  .  ,  1*0  .523  I 

(symmetric)  ^  0  J  oils 

X1  =  2.857,  e1  =  [.464,  .457,  .470,  .421,  .42 1  ]’ 

=  .809,  e2  =  [.240,  .509,  .260,  -.526,  -.582]' 

■*» 

Xj  =  .540,  e3  =  [-.621,  .178,  .335,  .541,  -.435]’ 

Xj  =  .452,  e4  =  [.387,  .206,  -.662,  .472,  -.382]’ 

A,.  =  .343,  ec  =  [-.451,  .676,  -.400,  -.  176,  .385]' 
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It  is  blear  we  have  two  distinct  groups  of  stocks  with  roughly  the  same 
within  group  correlation  (r  =  .55)  and  (very)  roughly  the  same  between  group 
correlation  (r  =  .40).  The  first  two  sample  principal  components,  Y^  =  e^'Z  and 

«v 

Y2  =  e2'Z  explain  a  proportion  (2.857  +  .809)/5  =  .73  of  the  total  variance  of 
the  standardized  returns.  Clearly,  the  first  component  represents  a  "suin"  of 
the  five  stocks  and  the  second  component  represents  a  "contrast"  between  the 
chemicals  and  oils. 

Assuming  a  k=2  factor  model  of  the  form  (19)  for  the  stock  returns  and 
estimating  the  loading  matrix,  L,  with  L  =  [/x^e^,  we  might  identify  the 

*w. 

first  factor  as  a  "market  factor"  and  the  second  factor  as  an  "industry"  factor 
(see,  also,  Johnson  and  Wichern,  1982,  chapter  9).  All  stocks  have  large  and 
nearly  equal  loadings  on  the  first  factor.  The  second  factor  is  a  bipolar  fac¬ 
tor.  The  chemical  stocks  have  positive  loadings  and  the  oil  stocks  have  negative 
loadings  on  this  factor. 

Comparing  the  correlation  matrix,  R,  above  to  the  correlation  matrix  p  dis- 

•w 

cussed  in  the  numerical  example  of  section  2,  we  see  that  there  is  some 
(descriptive)  evidence  that  the  correlation  patterns  among  stocks  are  somewhat 
regular.  Specifically,  if  we  consider  two  groups  of  stocks  with  m=3,  n=2  and  a 
population  correlation  matrix  of  the  form  (2)  with  p  =  .6  (the  case  considered 
in  the  numerical  example  of  section  2),  the  first  two  (important)  eigenvalues 
and  normalized  eigenvectors  of  p  are 

X1  =  2.832;  e1  =  [.469,  .469,  .469,  .412,  .412]' 

>2  =  .969;  e2  =  [.336,  .336,  .336,  -.575,  -.575]' 

■S, 

Wiese  are  very  similar  to  the  first  two  eigenvalues  and  eigenvectors  of  the 
sample  correlation  matrix  R.  Moreover,  the  three  smallest  eigenvalues  of  R  are 
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roughly  equal  ( X  =  .445  for  X^»  X^  and  X^)  and  consistent  with  the  final  three 
eigenvalues  of  p  ( X^  =  X^  =  X,.  =  .400). 

*>w 

We  have  provided  some  evidence  for  the  linear  factor  model  for  stock 
returns  assumed  by  Ross.  In  addition  we  have  suggested  that  the  measures  of  sys¬ 
tematic  risk  (the  factor  loadings)  may  be  more  regular  than  initially  proposed. 

6.  Summary 

We  have  derived  general  expressions  for  the  eigenvalues  and  eigenvectors  of 
certain  patterned  correlation  matrices.  The  particular  patterns  examined  were 
suggested  by  observed  correlations  among  common  stock  returns.  We  have  discussed 
our  results  in  the  context  of  a  principal  components  or  factor  analysis  of  stock 
returns.  The  limited  evidence  presented  here  suggests  the  linear  factor  model 
for  stock  returns  suggested  by  Ross  as  part  of  the  Arbitrage  Pricing  Theory  of 
assets  is  not  unreasonable.  Moreover,  the  factor  loadings  (measures  of  the 
response  of  the  returns  to  the  systematic  risk  components  of  the  model)  may  be 
equal  for  distinct  groups  of  stocks.  This  implies  more  regularity  in  the  gener¬ 
ating  mechanism  for  stock  returns  than  previously  imagined.  However,  more  work 
is  needed  before  this  assertion  is  confirmed. 

Finally,  we  note  that  the  derivations  for  the  matrix  p  (see  (21))  presented 
in  section  4  can  be  extended  to  the  case  of  several  matrices  B  and  C  appearing 
in  each  row.  However,  the  numerical  calculations  involved  become  rather  cumber¬ 
some. 
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